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Abstract: We present a way to include impurities in AdS/CFT correspondence, in 
view of its application to condensed matter physics. Examples of these are the current 
impurity and spin impurity We calculate electric conductivity and spin susceptibility 
of holographic superconductors, with doping of density/spin impurities. 
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1. Introduction 

Holography or AdS/CFT correspondence [1, 2, 3] is an extremely useful tool, in the 
sense that it gives a new view of the hard-core problem of quantum field theory, 
strongly coupled limit of the theory, from the gravitational theory or string theory 
viewpoint. Recently there are many advances applying the holography to more real- 
world setup like nuclear physics, hadron physics, and condensed matter physics such 
as quantum hall effects and superconductivity. For the case of condensed matter 
application, in many of the holographic setting, even though its validity to real- 
world condensed matter system is not clear, there is a big hope that these totally 
new holographic perspectives will shed new light on long-standing problems of these 
fields in the end. And in fact, recent progress along these directions, especially 
condensed matter physics, is quite remarkable. For reviews of recent development 
along these fields, see for example [4, 5, 6, 7, 8, 9, 10]. 

Many of the transport coefficients in strongly coupled systems can be calculated 
by using holographic setting. In condensed matter physics, there is a crucial in- 
gredient, the effects of impurities. Adding impurities is very ubiquitous technique 
in condensed matter physics, it changes the basic nature of materials and changes 
the transport properties drastically. In other words, it changes the Hamiltonian 



of the system. In a very famous example of high T c cuprate superconductors, by 
compounding Sr or Ba to La-Cu-0 system to add the conductance charge carriers, 
the properties of the material drastically changes from a Mott insulator to a high 
T c superconductor. Another famous example is the Kondo effect, where magnetic 
impurities, i.e., spins, make the resistance increase logarithmically as we lower the 
temperature. As is seen by these drastic changes of the properties of the material, 
impurities can induce many peculiar effects to the material. Therefore, it is very nat- 
ural to ask, how the impurities give effects generically, in many of the holographic 
condensed matter settings which were studied before. In this paper we will take the 
first step for this direction, especially the effects of the impurities, which include 
magnetic impurities. 

Kondo effect is induced due to the quantum spin nature of magnetic impurities 
coupled to conducting electrons near their Fermi-surfaces. The spin operator of the 
impurity fermion t/Vl 1 can be described by 



In this paper, treating as an external input, we consider the effects of the impurity 
to the transport coefficients in holographic settings. Our impurity is introduced 
by hand. For simplicity, we consider only homogeneous impurities which induce 
homogeneous source current J M . This naturally induces a bulk current J M (r), which 
is dependent only on the AdS bulk radial coordinates r. 

The organization of this paper is as follows; In section 2, we first describe in detail 
our method to include impurities in the AdS/CFT correspondence. Then in section 

3, we use the holographic superconductor model [11] and introduce the impurity. We 
calculate the AC conductivity, and evaluate the change of it due to the impurity. We 
treat two kinds of impurities there: density impurity and spin impurity. In section 

4, we calculate the spin susceptibility in the same manner. In both cases, we clearly 
see the effect of the impurity. Section 5 is for the discussion. 



Impurities play an indispensable role in condensed matter physics, while its treat- 
ment in the AdS/CFT approach has not been fully established. There are several 
interesting works [12, 13] which introduced local impurities, but in general, impuri- 
ties can be distributed all over the material of concern. Here we present our generic 
method to include some impurities in a general setup in AdS/CFT correspondence 2 . 

1 This impurity fermion ip- lm should not be confused with conducting quasi-particle fermions, like 
electrons. 

2 See also [14], where the authors studied the momentum relaxation due to impurities. 




2. Impurities in AdS/CFT 
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2.1 How to include the impurities 

Without the impurities, the conductivity (or the current correlator) is treated in the 
AdS/CFT correspondence as follows. Let us consider a strongly coupled fermion 
i/j in 1 + 2 dimensions. The fermion number current is given by = ip'j^if). We 
weakly gauge this U(l) symmetry, which is nothing but the electromagnetism, to 
have the coupling J (ftxA^O^ . To derive the current correlator, in the gravity side, 
we upgrade the "source" A^x^) for the operator to a bulk field Am{x^ ^ r) where 
r is the bulk coordinate along the AdS radial direction. According to the AdS / CFT 
dictionary, we demand that the bulk field approaches to the value of the boundary 
source, lim^oo A M=fM = A^. We solve the equation in the bulk, and substitute it to 
the bulk action, to obtain the classical bulk partition function which is a function 
of the boundary source A^. In the AdS/CFT [2, 3], this is equal to the boundary 
partition function with the source A^ of the source term J d^xA^O^. 

Now, let us add the impurities. Suppose the impurity is giving the electric 
density and the electric currents only, then the impurity coupling at the boundary 
theory is the minimal coupling J d^xA^O^ , where the impurity operator would be 
given by the impurity fermion field ipi m (x^) as = < 0i m 7 , Vi m . Let us treat the 

impurities as classical objects. Since we can distribute the impurity in an arbitrary 
manner in real-world experiments, we can take the vacuum expectation value of the 
current (O-. l mp (x fJ ')) to be an arbitrary function. That is, our impurity coupling is 



In going to the gravity side of the AdS/CFT correspondence, we upgrade this 
coupling (2.1) to the bulk coupling, 



We propose that this is a generic effect due to the impurities in AdS/CFT correspon- 
dence. If one specifies the dynamics of the impurity, the configuration J (a; M ,r) is 
determined. Basically, the radial dependence of the bulk source field J M {x^,r) rep- 
resents how the impurity responds in different energy scales. Since the behavior of 
the impurities at different energy scale can be taken arbitrary as an external input, 
we take J M as an input source to the bulk gravitational action. 

To gain more insight on the reason of this bulk source coupling, let us consider the 
impurity in the following manner. First, generic impurities are things put by hand. 
So we may allow such arbitrariness. Second, impurities can be heavy compounds 
and then their motion is not taken into account for transport coefficients. So the 
effect of the impurities are fixed in the bulk, and affects directly the currents of the 
conducting electrons (which are gauge bosons in the bulk). This allows us to treat 




(2.1) 




(2.2) 
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J as a classical input. In section 2.3, we provide explicit models to derive the 
coupling (2.2), as an example. 

In this way we consider the effect of the impurity in the form of the function 
J M (x /1 ,r) which we give by hand. 

2.2 Spin impurities 

In most interesting cases, impurities couple to the electromagnetic fields (and to the 
conducting electrons and quasiparticles in strongly correlated systems) through a 
spin magnetic coupling. In the previous subsection we have treated a covariant and 
canonical coupling between the impurity fermion and the electromagnetic field. 
Let us consider the spin magnetic interaction, which is written as 

J d 3 x BiWi&iaf)), (2.3) 

where 5 l (x M ) is the spin operator made from the impurity fermion field and 
Bi(x) is the magnetic field. 3 In terms of the field strength, this coupling is written 

as 

J d 3 x l -e ljk F jk {x)(S\x»)) = - J d 3 x A fc (x)(e^SV)>, (2-4) 

which is of the form (2.1). So, introducing the spin magnetic coupling for the impurity 
is included in the scheme of the previous subsection. 

In this paper, in order to show explicit examples, we consider a certain specific 
form of the impurity source J (x M ,r) in the bulk (which we give by hand), namely, 

J M = c/r n , (2.5) 

and compute the AC conductivity in the x direction. Obviously one can generalize the 
functional form of the impurity and also the direction of the currents. For numerical 
calculations, we specifically treat the case of n = 6 as explicit examples. We have 
to require that the back reaction of the source J may not change the asymptotic 
geometry, so we choose the source J which quickly decays near the boundary r = oo. 
The radial r dependence in the geometry corresponds to the energy dependence of 
the impurity source. 4 

According to (2.4), the source J y is related to the spin roughly as J y ~ d x S z . 
So, if we define a bulk version of the source spin, it is (for our choice of J as (2.5)) 
S z ~ x/r 6 . For numerical purpose, in this paper we consider only this case in section 
3 and 4. 

3 Here for simplicity the material is supposed to have 2 spatial dimensions, but one can generalize 
this to 3 spatial dimensions. Note that even in 2 spatial dimensions, on materials one can introduce 
3-dimcnsional spins. 

4 If one regards the energy as a hypothetical temperature of the impurity, the asymptotic decay 
of the source J means that near the boundary, where the energy (temperature) scale is large, the 
effect of the impurity is lower. 
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2.3 Impurity models in holography 

In section 2.1, we have argued that the way we introduce the impurity is just the 
introduction of the source term in the bulk, (2.2). Here, for concreteness, we shall 
provide explicit AdS/CFT models in which the impurity supplies the bulk source 
term. 

Let us consider a theory in 1+2 dimensions, with a strongly correlated fermion ijj 
and a strongly correlated massive impurity fermion ip- im . According to the gauge/gravity 
correspondence, we consider operators of our concern, 

o» = fty^, OL = fey'Vim- (2-6) 

These are just fermion number current for each fermion. 5 The AdS/CFT correspon- 
dence requires bulk fields which corresponds to these fields, 

0" Ap, 0f m ^ Bp. (2.7) 

Now, because these two currents are conserved, we write a gauge-invariant action for 
these two U(l) gauge fields: 

L = -\f^ - l -H^ + cH, u F^, (2.8) 

where 

Ffiv = d fl A u - d v A„ Hp, = d, L B u - d u B^. (2.9) 

The last term in the Lagrangian is one example of a coupling between the conducting 
electron and the impurity fermion. 6 In general, we can allow arbitrary couplings 
which are not prohibited by the symmetries, and this is the lowest order term in the 
derivative expansion. 

If impurities are heavy, we can generally fix the dynamics of the impurity-induced 
gauge potential £> M . So we can suppose that we may consider a background config- 
uration for the bulk field B^. Once we allow it, substituting this B yi configuration 
back into the action (2.8) and make an integration by parts, then we arrive at the 
bulk source coupling (2.2), A^J 11 in the bulk. 7 

5 For example, one can imagine an SU(N) gauge theory with two kinds of fermions ip and ipi m , 
both of which are in the fundamental representation of the SU(N). Then the model looks a QCD 
with two flavors. If we suppose that the fermion tp is light while the impurity fermion ipi m is 
superheavy, then the model suffices our purpose. In this paper, we would like to consider more 
generic situation in holography. 

6 The coupling is similar to the one for realizing a vector meson dominance in hidden local 
symmetry models. We would like to thank K. Fujikawa for pointing this out to us. 

7 Strictly speaking, the electromagnetic U(l) is a linear combination of the first and the second 
U(l). However, since the impurity dynamics can be killed due to its heaviness, it is equivalent to 
consider only the as dynamical field to calculate the conductivity 
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The configuration of the source J in the bulk depends on the dynamics of the 
impurity. For example, the asymptotic behavior of the source J would be related to 
the form of the coupling and also the conformal dimension of the impurity operator. 
Furthermore, we may have a variety of the form of the couplings and also we may 
allow multiple kinds of the impurity operators with different conformal dimensions. 
This means that basically an arbitrary configuration is allowed for the bulk source 
J, depending on the models. Therefore in this paper we do not stack to a specific 
model, and we allow generic configuration of the bulk source J. The idea is along 
the standard picture that one can control the distribution of the impurity source by 
hand. 

We present one more model which would suffice for introducing a spin back- 
ground according to the idea written above. The impurity spin operator is 

SL = ^m7¥Yim, (2.10) 

in 1+3 dimensional notation for fermions. This can be easily seen if one decomposes 
the four-component spinor into 2+2 components (<p, x) 7 \ and explicitly write 
the gamma matrices in terms of the Pauli matrices: S l = tp<a i tp + x a% X- So we can 
consider the spatial component of the axial current as a spin density. Let us introduce 
a bulk vector field Bp which is dual to the operator S l , in the same manner. Then, 
since the axial current is not conserved for the massive fermions, the corresponding 
£7(1) gauge symmetry for Bp is broken in the bulk. Therefore we can write a model 
like 

L = --^F^ - l -H, v H^ + l -m 2 B,B» + cH, u F^. (2.11) 

Note that we have the mass term which breaks the U(l) gauge invariance for the 
field Bp, so Bp is a massive gauge field (Proca field). The model, with the coupling 
c, explicitly introduces the interaction between the conducting electrons and the 
impurity spins. The mass term can be replaced by a condensation of a charged 
scalar field in the bulk, as in the case of the holographic superconductors [11]. 

The spin is a pseudo- vector, as the current (2.10) is a pseudo- vector. In the 
model (2.11), the last term breaks the parity invariance. To avoid such a violation 
of the parity symmetry, one can write other kind of the coupling. For example, the 
following higher derivative coupling 

ce^Fp U H pa B x B x (2.12) 

preserves the parity symmetry. (Note that without the last factor B\B X the coupling 
is trivially zero through an integration by parts.) 

It is of course possible to proceed with these explicit models more, but it is not 
our purpose of this paper: we argue that in this way generic effects of the impurity 
can be taken into account as a bulk source term. In the rest of this paper, we just 
take a single configuration of the source J which is (2.5) with n = 6, as an example. 
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3. Conductivity 



In this section, we calculate the conductivity of a superconductor, in the presence of 
the impurity. As we studied in the previous section, the impurity can be introduced 
via the source term for the bulk gauge field. Once we include the source term, the 
equations of motion for the bulk fields are modified. 

In general, the introduction of the source term changes not only the gauge field 
configuration but also the background geometry. It is desirable to calculate the back- 
reaction to the geometry, but in this paper we shall not consider it, by taking the 
gravitational coupling constant to zero. In other words, we take a probe limit, as the 
gauge sector works as a probe for the geometry. The limit is useful in two senses: 
first, it can pick up intrinsic physics of the electromagnetic currents as they live in 
the gauge field sector, and second, it offers a simpler ground for calculations. 

First, we solve the background configuration of the gauge fields and the scalar 
fields, with the source term which represents the impurity. Then, we study the 
fluctuation to calculate the conductivity. Our calculation follows what has been 
originally proposed in [11] for the holographic superconductors. 

3.1 General strategy 

3.1.1 The background configuration in the bulk 

The Lagrangian in the bulk is that of a Maxwell field and a charged scalar field, 



S = d g 



(3-1) 



Note that we have added the last term, the source term, which represents the impu- 
rity. The scalar potential (2/L 2 )|\l/| 2 can be of a different form, but we just follow 
the popular example presented in [11]. 

The background metric is the Schwarzschild black brane, 8 

ds 2 = -f(r)dt 2 + -^-dr 2 + r 2 (dx 2 + dy 2 ), (3.2) 
f( r ) 

r 2 M 

f(r) = j- 2 --- (3.3) 
L z r 

The location of the horizon r = r + and the temperature T of the black brane is 



8 One may wonder if one can take an AdS-Rcissncr-N6drstrom black brane as a background 
metric, as anyway we are going to consider a nontrivial configuration of Aq(t). However, with the 
Rcissner-Nodrstrom background, in the probe limit k — > 0, the temporal component of the gauge 
field should go to infinity. Then it is difficult to make sense of the the fluctuation of the gauge field 
around such a divergent background. 
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In this background metric, we solve the equation of motion of (3.1), with a given 
source J(r). 

The equation of motion for the gauge field and the scalar field is 



1 



d^yf^Fn + ig vX (V*(d x - iA x )V - V(d x + iA x )**) = .F , (3.5) 



V~9 

,(d, - iAJ W=gg» v {d v - iA v )V] + = . (3.6) 



We need to solve these equations to obtain the background configuration of the gauge 
field and the scalar field. 

3.1.2 Calculation of the conductivity 

Given the background solutions of equations (3.5) and (3.6), we can now calculate 
the conductivity. For the calculation of the conductivity, say, along the x direction, 
we turn on a fluctuation of A x which is time-dependent with frequency u as ~ e~ tujt . 
The equation for this A x fluctuation is 

K + +(^- 2 -^-)a x = 0. (3.7) 



/ ' V/ 2 / 

Note that this is the same equation as what is given in [11], since the source term 
which we introduce affects only the background configuration. 9 The boundary con- 
dition to solve the fluctuation equation is the in-going boundary condition at the 
horizon. 

As explained in [11], the conductivity a(u) is calculated by the ratio of the 
asymptotic coefficients of the normalizable and the non-normalizable modes of A x . 
Expanding 

A x = Af> + — ^- H (r ~ oo) (3.8) 

r 

the conductivity is given by 

a(co) = f^. (3.9) 

uAr 

In the following, we treat two examples: (i) the source is J*, meaning that the 
impurity affects the electric density, and (ii) the source is J y , meaning a particular 
kind of the impurity spin density. 

9 If one introduces a source term which is not linear in A x but quadratic in A x , it will affect 
directly to the fluctuation equation. The choice of the form of the term with the impurity is 
dependent on what kind of interaction one wants to introduce in the boundary theory. 
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Figure 1: The plot of the background configuration of $(r). Thick line is with the 
impurity J* = 0.1/r 6 , and the thin line is without the impurity. 

3.2 The case of impurity density 

First, we consider the example with J*(r). Obviously, for the background configura- 
tion of the gauge field, we need only the temporal component, 

A t = $(r). (3.10) 

We will take a gauge where \I/ takes the real value. 

The equations of motion for $(r) and \P(r) are easily derived as 

2 2^ 2 

$" + -$'-—$ = J\ (3.11) 
r f 

(V 2\ $ 2 2 

+ {j + r) *'+p* + j2]-* = °- ( 3 - 12 ) 

The boundary conditions we need to care about are $(r = r + ) = at the horizon, 
and at the asymptotic boundary r ~ oo, 

$(r) £ + -.. , (3.13) 

r 

^(r) ~ — + — + ... . (3.14) 

As chosen in [11], we put either or \E , ^ 2 ' ) to vanish. For simplicity in this paper, 
we put \l/( 2 ) = 0. Then the input variable is only the charge-career density p (and the 
temperature T of the background geometry) . The chemical potential \i is determined 
by p. 

For a given form of the source J*, it is easy to perform a numerical calculation for 
$(r) and \&(r). The plot of the resultant configuration, with and without the source 
term, is shown in Fig. 1. For the numerical parameters, we have chosen J* = 0.1/r 6 , 
and p — 1, L — 1, and M = 0.2. The temperature is T ~ 0.14 which is in the 
superconducting phase [11]. 

Next, we turn to solve the fluctuation equation (3.7) to obtain the conductivity. 
It is straightforward to solve the fluctuation equation numerically, and the result for 
the real and the imaginary parts of the conductivity is shown in Fig. 2. 
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Figure 2: The plot of the conductivity a as a function of the frequency uj. Thick lines: 
conductivity with the impurity J* = 0.1/r 6 . Thin lines with dots: conductivity with no 
impurity. Left: Real part of a. Right: Imaginary part of a. 

Interestingly, the real part of the conductivity increases due to the impurity. 
A possible interpretation of this effect is as follows: since we have introduced the 
impurity as a background density J*, we may think of it as an increase of an "effective 
mass" of the background geometry. Once the mass of the black brane increases, the 
effective temperature increases. This generally increases the conductivity, which is 
consistent with our numerical result. 

3.3 The case of impurity spin density 

We provide another example, which is J y (r). As explained earlier, turning on the J y 
component represents a background spin configuration caused by the impurity. 

With this source term in the bulk, the A y component is excited as a background. 
The equations of motion are 



The boundary condition for the A y is A y ~ 0(l/r) at r ~ oo, as the background 
source J y should induce a normalizable mode of A y . 

The numerical calculation of the background gauge and scalar configuration, and 
also of the fluctuation (3.7), is straightforward, and the results are shown in Fig. 3 
and Fig. 4. Interestingly, the result looks quite similar to what has been obtained for 
the case of J*. The interpretation would be the same as that for J': the background 
impurity may cause an effect of an increase of the "effective temperature" . 

Note that if we measure the conductivity along the y direction, the result is 
different. This is because the fluctuation of A y couples 5^, as the background A y is 
nonzero. In addition, once we work without the probe limit, there should appear a 
difference between the conductivity along the x and the y directions. 




(3.15) 



(3.16) 



(3.17) 
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Figure 3: The plot of the gauge configuration (<$, Left) and the scalar configuration 
Right). Thick lines are with the impurity while the thin lines are without the impurity. 
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Figure 4: The plot of the conductivity a as a function of the frequency u, for the spin 
impurity J y (r). Thick lines: conductivity with the impurity J y = 0.1/r 6 . Thin lines with 
dots: conductivity with no impurity. Left: Real part of a. Right: Imaginary part of a. 



4. Spin susceptibility 

In this section, we calculate the effect of the impurity on the spin susceptibility. 
As we have described in section 2, the bulk source J v we have introduced may be 
considered as a spin impurity. So it should have some effect on the spin susceptibility 
in the measurement. 

The spin susceptibility can be measured by a correlator of magnetic fields. Its 
treatment in the context of the AdS/CFT correspondence was discussed in, for exam- 
ple, [15]. The magnetic field fluctuation is a kind of the gauge field fluctuation which 
we considered in the previous section for computing the AC conductivity. The only 
difference is the fact that we need a spatial modulation for to induce the typical 
magnitude of the fluctuation magnetic field: Concretely, we consider the fluctuation 
of the form 

5A X = e- iujt e igy a(r), (4.1) 

where q is the momentum along the y direction. This is equivalent to consider the 
fluctuation of the magnetic field, 10 

SB 3 = -d y 5A x = -iqe- lu)t e iqy a{r). (4.2) 
10 We may consider a different orientation for the spatial modulation. For example, even for the 
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Figure 5: Rigid thick line is the spin susceptibility with the impurity, while the thin line is 
for no impurity. We have chosen the impurity J y = 0.2/r 6 as before, and took the spatial 
modulation q = 1 for simplicity. The real part of the spin susceptibility increases for the 
middle region of the frequency uj under the doping of the impurity. 



The ratio of the normalizable and the non-normalizable modes gives the spin sus- 
ceptibility, so, the procedure to obtain the spin susceptibility is the same as that for 
the conductivity. 

The fluctuation should satisfy the on-shell equation. Substituting the fluctuation 
(4.2) and the background geometry to the on-shell equation, we obtain the following 
equation 

i n yy \]> 2 

TT^rJr {V=gg rr g xx d r a(r)) - q 2 ^a(r) - 2—a(r) = u 2 a{r) . (4.3) 



y—ggg g g 

The numerical calculation of the spin susceptibility is straightforward. The result 
is shown in Fig. 5. We observe that the spin susceptibility Xs increases due to the 
impurity doping. 



5. Discussions 



In this paper, we proposed a generic method to include impurities in the context of 
AdS/CFT in application to condensed matter physics. We have calculated several 
examples for the transport coefficients in a superconducting phase: AC conductivities 
and spin susceptibility. The effects are clearly seen. 

The essence is to introduce a source term in the bulk side of the AdS/CFT 
correspondence which sources the bulk gauge field. The bulk source term may have 
an arbitrary profile, which corresponds to the arbitrariness of how we introduce the 
impurities to the system. In section 2, we have given explicit models which can 
determine the profile of the bulk source term from the property of the impurity, 
but they are just examples: in principle, we may allow arbitrary profile for the bulk 
source term in general. 

same spin direction B 3 , we may take 8A y = e~ iUt e tgx a(r), which gives the magnetic field fluctuation 
5B 3 = d x SA y = iqe- iu]t e lqx a(r). 
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For the example which we analyzed in this paper, the impurity in the holographic 
superconductors, we observed that the AC conductivity and the spin susceptibility 
increase due to the impurity The former increases at the lower frequency, while 
the latter seems to be sensitive to the typical spatial momentum of the magnetic 
field. Since the impurity is introduced as a small perturbation, it is natural that the 
impurity affects only the low frequency region of the AC conductivity. The increase 
may be interpreted as an increase of an "effective temperature" of the system. We 
have done the calculation with different bulk profile for the bulk source field, and 
obtained a similar behavior. It could be that the increase at the low frequency region 
which we obtained is a generic feature, but we do not have a clear understanding for 
this. 

We have worked in a probe limit where the effects of the gravity is ignored. 
In that case, the background geometry unperturbed with the impurity source term 
would be a Reissner-Nodrstrom black brane with a scalar hair [16]. There are two 
effects by the gravity effects. One is the back-reaction to the background geometry 
due to the impurity. The other is the gravity fluctuation on the background, which 
couples to the gauge field fluctuation for the conductivity calculation. 

One can also introduce other fields, like dilaton and axion, in the bulk action 
to simulate quantum Hall systems as in [17, 18]. If we have running dilaton, the 
geometry can be replaced by Lifshitz type of geometry or others [19, 20, 21, 22] at 
IR where the geometry shows different symmetries, such as non-relativistic one. In 
this way, on the gravity side, we have varieties of generalizations. In the holographic 
superconductors, some impurity of course would change the critical temperatures, 
which would be interesting, and that could be calculated when one includes the 
coupling to the gravity. Various transport coefficients are affected by the impurities, 
while the universality of the critical exponents against perturbations by the impurities 
would be an interesting question. 

The impurity we considered in this paper is "classical" impurities, while in re- 
ality impurity may be caused by quantum impurity fields. If one allow quantum 
impurities, nontrivial dynamics of the impurities would determine the total effect to 
the conducting quasiparticles. For example, one-loop effect of the impurity fermions 
would become the leading order effect [23, 24]. 

The treatment of the bulk source field in this paper is in the Maxwell sector, so 
the effect of the impurity looks just a linear order in J. However, the coupling to the 
scalar field is nonlinear. Or, once one couples it to a dynamical gravity, higher order 
term in J will show up, due to the nonlinearity of the gravity. There are interesting 
effects which appear only at higher order in the impurity perturbation, such as the 
renowned Kondo effect. As we pointed out at the end of section 3, there should be 
anisotropy effects for the conductivity due to the non-linearity. It is interesting to 
investigate these furthermore. The exact treatment of the bulk source field in the 
dynamical gravity model is intimately related to the dynamics of the impurity, so 
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solving models beyond the probe limit is of importance. 
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